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(54) Method for transmission line impulse response equalisation and a device to perform this 
method 



(57) A set of parameters (W) is calculated for an 
equalizer which equalizes a transmission line impulse 
response (LIR) in such away that the equalized impulse 
response (EIR) approximates a target impulse 
response (TIR) with predetermined length (L). 

An error function (E) is minimized via eigenvalue- 
eigenvector calculation of a channel dependent matrix 
(O). This channel dependent matrix (O) is build up from 
a signal (X), a distorted signal (Y), and length (L) and 
delay (D) of the target impulse response (TIR). The 



error function (E) comprises a first contribution which 
represents the difference between the equalized 
impulse response (EIR) and the target impulse 
response (TIR) and a second contribution which repre- 
sents the energy (E DC ) that is transmitted in unused fre- 
quency bands. The eigenvector corresponding to the 
minimum eigenvalue of the channel dependent matrix 
(O) represents the set of equalizer parameters (W). 
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Description 



The present invention relates to a method for calculating a set of equalizer parameters as described in the pream- 
ble of claim 1 , a device for calculating such a set of equalizer parameters as described in the preamble of claim 9, and 

5 a demodulator including such a device as described in the preamble of claim 10. 

Such a method and equipment to perform such a method are already known in the art, e.g. from the article 'Adap- 
tive Channel Memory Truncation for Maximum Likelihood Sequence Estimation', written by D.D. Falconer and F.R. 
Magee, Jr., and published in The Bell System technical Journal, Vol: 52, No. 9, November 1973, from page 1541 to 
page 1562. As described in the second paragraph of the introduction on page 1541, an equalizer or linear prefilter 

10 which forces the overall impulse response of the channel/prefilter combination to approximate a desired truncated 
impulse response is presented therein. The overall impulse response of the channel/prefilter combination is equivalent 
to the equalized impulse response mentioned in the preambles of claims 1, 9 and 10, whilst the desired truncated 
impulse response in the above mentioned article is called a target impulse response with predetermined length in the 
present application. In the known method, the desired impulse response is obtained by calculating an eigenvector cor- 

15 responding to the minimum eigenvalue of a certain channel dependent matrix. This channel dependent matrix and 
more specifically the minimum eigenvalue thereof is representative for a mean square error between the output of the 
prefilter and the output of the desired prefilter, the desired prefilter having an impulse response equal to the just men- 
tioned desired truncated impulse response. When minimizing the mean square error, the desired impulse response is 
optimized. As is proven in section III of the above cited article, the known method then adaptively or iteratively calcu- 

20 lates the linear prefilter parameters from the desired truncated impulse response. Thereto, in successive iterations, the 
Viterbi algorithm is applied to the prefilter output, the result thereof is applied to the desired prefilter and the output of 
this desired prefilter is subtracted from the prefilter output in such a way that an error sequence for adjustment of the 
desired impulse response and prefilter parameter set is obtained. The receiver structure allowing to perform the known 
method is drawn in Fig. 3 on page 1545 of the above mentioned article. Although the known technique or derivatives 

25 thereof wherein e.g. the Viterbi algorithm is not used exhibit very low mean square error values and as a consequence 
a theoretically negligible intersymbol interference, its applicability in ADSL (Asymmetric Digital Subscriber Line) is 
doubtful. The main reason is linked to the finite precision limitations of the hardware. Indeed, the equalizer transfer func- 
tion obtained when applying the known technique in ADSL, will boost the unused frequency band, i.e. the frequency 
band below 25 kHz, and will attenuate the tones in the ADSL passband, which approximately cover the frequency band 

30 from 25 kHz to 138 kHz. A first consequence is the decrease of numerical stability. It is even possible that the attenua- 
tion is so severe that the first stages of an FFT demodulator as is used in ADSL receivers rounds off the frequency com- 
ponents to zero such that the frequency domain equalizer of the receiver is unable to boost these tones again. An 
additional consequence is that, by implementing the known method in ADSL, the attenuation is very severe at the filter 
boundaries between the unused and used frequency bands. This is a result of the very sharp filter edges there. The 

35 pilot tone will very likely be located close to the region where the transmission channel exhibits a relatively low attenu- 
ation. This is mostly located near the lower downstream filter edge. Synchronization may therefore be subjected to an 
unexpected decrease of the signal to noise ratio of the pilot tone 

An object of the present invention is to provide a method and equipment for calculating a set of equalizer parame- 
ters of the above known type but wherein numerical instability due to amplification of unused frequency bands and 

40 attenuation of used frequency bands is optimized, and wherein the above mentioned synchronization problem is 
avoided. 

According to the invention, this object is achieved in the method, device and demodulator described in claims 1, 9 
and 10 respectively. 

Indeed, by adding to the difference between the equalized and the target impulse response, a second contribution, 
45 an additional degree of freedom which allows to shape the equalizer impulse response is provided. By adding a new 
error source to the mean square error criterion, the present invention prevents boosting unused frequency bands 
instead of used frequency bands. The higher the energy located in the unused frequency band, the higher the additional 
error. As a result, the present method searches for equalizer parameters which define an equalizer impulse response 
that does not boost unused frequency bands. From the above explanation, it is clear that in the present invention, the 
so shortcomings with respect to numerical stability and synchronization are overcome. 

A further feature of the present invention is that the first contribution to the channel dependent matrix, representing 
the distance between equalized and target impulse response, has the form described by claim 2. 

Indeed, as will be proven later on in the description, modeling the transmission line in a certain predetermined way, 
allows to build up the first part of the channel dependent matrix O as given by the mathematical expression in claim 2. 
55 In different implementations of the present method, step d of generating a channel dependent matrix O may be sim- 
plified. Therefore, assumptions as described e.g. in claims 3 or 4 have to be made. Assuming, as described in claim 3, 
that the transmission line is not plagued by noise, or assuming, as described in claim 4, that the spectral density of the 
input signal applied to the transmission line by the transmitter has a flat shape, allows to simplify the mathematical 
expression of the channel dependent matrix O. and evidently reduces the number of calculations necessary in step d 
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to build up this matrix. For further details, reference is made to the description later on in this application. 

An additional characteristic feature of the present invention is that, in implementations thereof wherein the trans- 
mission line is not supposed to be noiseless, the noise is measured in a step included in the present method before the 
channel dependent matrix O is build up. This feature is described in claim 5. 

s A further characteristic feature of the present method is that, still in implementations thereof wherein the transmis- 

sion line is not supposed to be noiseless, the noise is taken into account in a way described in claim 6. Indeed, by mod- 
elling noise measured on the transmission line via an additive noise source, cascade connected to a noiseless 
transmission line, elements of the autocorrelation matrix RYY are modified as described in claim 6. Later on, it will be 
seen how the expression of RYY given in claim 6 is mathematically derived. 

io Yet a further characteristic feature of the present method is that, as described in claim 7, the second integer number 
of sample periods can be expressed as a function of the amplitude and the sample index number of the transmission 
line impulse response sample with highest amplitude. This feature corresponds to the intuitive feeling that the shorter 
target impulse response will approximate the transmission line impulse response in an optimized way only if this shorter 
target impulse response is delayed in such a way that the location in time of its peak values approximate the peak val- 

is ues of the transmission line impulse response. 

In a further implementation of the method according to the present invention, described in claim 8, steps c, d and 
e are executed for a sequence of acceptable delay values. In this way, for each delay value within an interval of accept- 
able delay values, a corresponding channel dependent matrix having its specific minimum eigenvalue is found. It will 
again be proven later on that this minimum eigenvalue is a measure for the error function that has to be minimized, and 

20 therefore, in this further embodiment, only the channel dependent matrix with smallest minimum eigenvalue is main- 
tained to calculate therefrom in steps f and g the set of equalizer parameters. 

The above mentioned and other objects and features of the invention will become more apparent and the invention 
itself will be best understood by referring to the following description of an embodiment taken in conjunction with the 
accompanying drawings wherein: 

25 

Fig. 1 represents a block scheme of an embodiment of a demodulator according to the present invention; and 
Fig. 2 represents a digital model for the cascade connection of a transmission line and an embodiment of an equal- 
izer whose time domain equalizer parameters are calculated according to the present method. 

30 In a theoretically ideal multicarrier system, i.e. in a system with infinitely large data symbols, the impulse response 
length of the transmission line is negligible compared to the length of data symbols transmitted over this transmission 
line. In known multicarrier systems however, acceptable latency or maximum acceptable delay between input and out- 
put of data symbols transmitted over a transmission line limits the length of these symbols. As is specified in section 
6.2.11 on page 22 of the draft American National Standard for Telecommunications on ADSL (Asymmetric Digital Sub- 

35 scriber Line), published by ANSI (American National Standards Institute) in April 1994, the length of a discrete multi 
tone (DMT) symbol is limited. In ADSL applications, the transmission line impulse response length is not negligible 
compared to the symbol length and intersymbol interference will occur. As a consequence, successively transmitted 
symbols are spread out and interfere as a result of which the transmission is distorted. 

In general, two different techniques are known to compensate for such intersymbol interference. In a first technique, 

40 called cyclic prefix extension, data symbols are extended at the transmitters. Gaps are created at the edges of the data 
symbols e.g. by inserting spare time between successive symbols. In these gaps intersymbol interference can occur 
without distorting the transmission. The transmission efficiency or throughput efficiency however, decreases with 
increasing lengths of the gaps or cyclic extensions. In the second intersymbol interference compensating technique, the 
transmission line impulse response is equalized by an equalizer, which is a digital filter at the receiver's side. The objec- 

45 tive of the equalizer is to compress the transmission line impulse response to a shorter overall impulse response. Full 
equalization of the transmission line impulse response however requires a large amount of computation. The draft 
ADSL Standard already mentioned above specifies the combination of both techniques. Equalizing the transmission 
line impulse response renders the transmission line impulse response length short enough to enable intersymbol inter- 
ference compensation via cyclic prefix extension with tolerable prefix length, and vice versa: compensation via cyclic 

so prefix extension allows the equalizer to be of acceptable computation complexity. 

In the following paragraphs, an algorithm will be derived which enables calculation of equalizer parameters in a fast 
and stable way. Before entering the mathematical approach of the problem, the ADSL (Asymmetric Digital Subscriber 
Line) environment, i.e. the demodulating part, wherein such an equalizer is used, will be described. 

The demodulator DEMOD of Fig. 1 includes an analog to digital converter ADC, a time domain equalizer TEQ, a 

55 serial to parallel converter SPC, a fast fourier transform processing unit FFT, a frequency domain equalizer FEQ, and a 
demapper DM which are cascade connected between an input Dl and an output DO of the demodulator DEMOD. 

Distorted signals at the output of the transmission line, connected to the demodulator input Dl but not shown in Fig. 
1 , enter the demodulator DEMOD via this input Dl. Such a distorted signal is sampled by the analog to digital converter 
ADC to be converted into a time domain digital sequence which is applied to the equalizer TEQ. The time domain equal- 
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izer TEQ is a digital filter whose set of parameters is calculated based on channel dependent information in accordance 
to the method of the present invention. In this way, by digitally filtering the time domain digital sequence, the overall 
impulse response length is shortened. The so obtained equalized time domain digital sequence is transformed into an 
equalized parallel time domain digital sequence by the serial to parallel converter SPC. The fast fourier transform 

5 processing unit FFT converts the equalized parallel time domain digital sequence into an equalized frequency domain 
parallel digital sequence and applies this sequence to the frequency domain equalizer FEQ which restorates amplitude 
and phase of the transmitted vectors, and to the demapper DM. This equalized frequency domain parallel digital 
sequence is a sequence of real data which represents a set of modulated carriers. These carriers have equidistant fre- 
quencies, carry different amounts of data bits which are modulated thereon via different modulation techniques. To 

w determine the amounts of data bits that are allocated to different carriers and to determine the modulation techniques 
used to modulate different carriers, a mapper included in the transmitter has to communicate with the demapper DM in 
the receiver. This mapper e.g. assigns 2 data bits to the first carrier and decides to modulate these 2 bits thereon via 4 
QAM modulation, assigns 4 data bits to the second carrier and decides that these 4 data bits are modulated thereon 
via 16 QAM modulation, and so on. The demapper DM included in the demodulator DEMOD at the receiver's side as 

15 aware of the modulation techniques used for the different carriers and of the amounts of bits allocated to each of the 
carriers. The task of the demapper DM is thus to demodulate each carrier to obtain therefrom the exact amount of data 
bits modulated thereon. The demapper DM then constitutes an output sequence and supplies this output sequence to 
the demodulator output DO. 

The block scheme of the Discrete Multi Tone demodulator DEMOD in Fig. 1 will not be described in further detail 

20 since such a detailed description is of no importance for the present invention. Further details with respect to ADSL 
requirements are described in the already mentioned draft ANSI Standard on ADSL whilst specific implementations of 
Discrete Multi Tone modulators and demodulators are found in the articles 'A multicarrier E1-HDSL Transceiver System 
with Coded Modulation', written by Peter S. Chow, Naofal Al-Dhahir, John M. Cioffi and John A. C. Bingham and pub- 
lished in the issue Nr. 3 May/June 1993 of the Journal of European Transactions on Telecommunications and Related 

25 Technologies (ETT), pages 257-266, and 'Performance Evaluation of a Multichannel Transceiver System for ADSL and 
VHDSL Services' from PeterS. Chow etal., published in the issue Nr. 6 August 1991 of the Journal of European Trans- 
actions on Telecommunications and Related Technologies (ETT), pages 909-919. 

In the following paragraphs, the method according to the present invention will be mathematically derived. Starting 
from a time domain model, a mathematical approach of the problem will lead to an expression of an error function and 

30 a channel dependent matrix which form the basis of the present method. 

The time domain model is drawn in Fig. 2. Therein, the transmission line is represented by the cascade connection 
of a noiseless transmission line with impulse response LIR and an additive noise source wherein noise N is added via 
a summator SM to the output signal of the noiseless transmission line. This output signal enters the summator SM via 
a first input 11 , whilst the noise N enters the summator SM via a second input i2 thereof. An input signal X applied to the 

35 transmission line appears at the output thereof as a distorted signal Y, distorted due to the transmission line impulse 
response LIR and due to the presence of the noise N. To compensate for these disturbances, a digital filter filters the 
distorted output signal Y This digital filter is included in an equalizer characterized by a set of parameters W. The trans- 
mission line impulse response LIR and the equalizer parameters W constitute an equalized impulse response EIR or 
overall impulse response. 

40 As already mentioned when describing the problem of intersymbol interference above, the global impulse 
response, i.e. the equalized impulse response EIR, should be short in length to enable avoiding intersymbol interfer- 
ence. Therefore, the equalized impulse response EIR in the digital model of Fig. 2 is compared to a target impulse 
response TIR, not further specified, but with a predetermined length L representing the guardband length or cyclic pre- 
fix length. This target impulse response TIR, delayed over a delay D has to approximate the equalized impulse 

45 response EIR to guarantee negligible intersymbol interference. Consequently, the mismatch between the equalized 
impulse response EIR and the target impulse response TIR has to be minimized. As will be seen in the calculations later 
on, an error E between the output Z1 of the equalizer and a fictive output 22 of the target impulse response TIR consti- 
tutes the first contribution to an error that has to be minimized. These outputs Z1 , and 22 therefore in Fig. 2 are applied 
to a non-inverting input Nl and an inverting input II of a subtractor SB respectively. 

so Following definitions are used in the mathematical approach later on in the description: 
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20 



25 



A(z) = 2> z " 



M 



H = 



h 1 

- h M-1_ 



the Z-transform A(z) of a digital sequence a, 
whose index i covers a range from -oo to +oc; 
the length of the transmission line impulse 
response LIR; 

the number of taps in the digital filter of the 

equalizer or, in other words, the number of 

time domain equalizer parameters W; 

the length of the target impulse response 

TIR; 

the delay applied to the target impulse 
response TIR; 

the line impulse response vector; 



30 



35 



40 



45 



50 



55 



BNSDOCID: <EP Q768778A1 J_> 



EP 0 768 778 A1 



w = 



10 



75 



20 



25 



30 



B = 



w 0 

Wi 

. W T-1 

bo 
b» 



x k-1 



|_ X k-IWMj 



k-O 
X k-D-1 

*k-D-L+1. 

y k .i 



35 



.y<c.T.ij 

Wopt 
Bopt 



40 



the equalizer time domain vector; 



the target, impulse response vector; 



a signal vector related to the line impulse 
response; 

a signal vector related to the target impulse 
response; 



a distorted signal vector related to the 
equalizer time domain; 



an optimal equalizer time domain vector; 

an optimal target impulse response vector; 

and 

a channel dependent matrix. 



45 



SO 



55 



When an input signal X is applied to the transmission line, a distorted signal Y appears at the output thereof. The 
k' th sample, k being an integer number, of the distorted signal Y does not only depend upon the k' th sample of the 
input signal X but also upon M-1 previous samples of this input signal X. This is due to the transmission line impulse 
response LIR which has length M. The transmission line is not perfect as a result of which input samples transmitted in 
the past still have an influence on actual output samples. The k' th output sample y k of the transmission line can be 
expressed as follows: 



M-1 



+ n. 



0) 



i=o 



Herein, i represents an index and the influence of noise on the transmission line is concentrated in the term n k , this 
being a result of the digital model of Fig. 2. Using the above defined vectors, the expression in (1) can be simplified as 
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follows: 



Yk = H T X hk + n k (2) 

5 

Herein, H T represents the transposed vector of H. 

The distorted output signal Y is applied to the equalizer to be digitally filtered thereby. Since the digital filter has T 
taps, the k' th sample z1 k of the equalized output signal 21 is given by the convolution 

10 T-1 

z1 k = £ w i-y k -i O) 

!=0 



is or, when using vector expression: 

z1 k =W T Y W)k (4) 



20 W T represents the transposed vector of W herein. 

When the input signal X further is applied to a fictive filter or fictive transmission line with an overall impulse 
response TIR whose length is restricted to L taps, a second, fictive, output signal 22 is obtained. This overall impulse 
response TIR is called the target impulse response and is defined by the target impulse response vector B. The shape 
of the equalized impulse response EIR which comprises the line impulse response LIR and the digital filter in the time 

25 domain equalizer W response will be determined in such a way that it approaches this target impulse response TIR. To 
allow the equalized impulse response EIR to approach the target impulse response TIR, a delay is applied to this target 
impulse response TIR. This delay should be considered as an additional parameter which has to be determined some- 
where in the method. The delayed target impulse response transforms the input signal X into the second output signal 
Z2, whose sample with index k is given by the convolution: 

30 

L-1 

z2 k = £b i .x k _ D _ i (5) 
i«0 



35 

or in matrix notation: 



22 k = B T X bk (6) 

40 

Therein, B T denotes the transposed vector of B. 

The target impulse response TIR, the time domain equalizer parameters W and the delay D now have to be deter- 
mined in such a way that the equalized impulse response EIR approaches the target impulse response TIR as close as 
possible. Therefore, output signals Z1 and Z2 are compared and the mismatch between both is minimized. This error 
45 is mathematically expressed by a sequence E whose k* th component e k is given by: sequence E whose kth component 
e k is given by: 

T-1 L-1 

e k = 21 k 22 k - I>i-y k-i - EMjmm (7) 

50 i=0 ioO 



or, when using the vector notation: 

55 



e k = W T .Y W(k -B T X bik (8) 



Thus, the error function that has to be minimized is equal to: 
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E{e 2 k } = E{(W T Y w>k -B T X bjk ) 2 } 



O) 



Therein, E denotes the mean square error between Z1 and Z2. This mean square error is not dependent upon k and 
s can thus be rewritten as follows: 



E{e 2 } = E{(W T Y w -B T X b ) 2 } 

io Since 

(W T . Y w -B T .X b ) 2 = W T .Y w .Y^.W-2W T .Y w .X b .B + B T .X b .xJ. B 

75 and by introducing the following correlation matrices: 

Ryy = E{Y W .Y*} 

Ryx = E{Y w .xJ} 

^xx = E{X b .xJ} 

(10) can be rewritten in the form: 

E{e 2 } = W T R W+B T R X x B " 2B T - r yx- w 



20 



25 



Herein: 



30 



35 



40 



Ryy = 



R»<(o) R„(1) 

RxxO) R»(o) 

R«(L-1) ... 

Rw(o) R W (D 

R„0) R w (o) 

Ryy(T-l) 



... R^L-1) 

RxxO) 
RxxO) Rxx(o) _ 

... R w (1) 

RyyCO Ryy(O) . 



45 



50 



Ryx(D) RJD + 1) 
Ryx(D-l) R„(0) 

R^D-T+1) 



Ryx(D+L-1) 



.. R^D + L-T) 



(10) 



(11) 



(12) 



(13) 



(14) 



55 



if it is assumed that the noise N is uncorrelated with the transmitted signal X. In (14): 

R xx (k)=E{x l+k .x,} 
Ryy(k)=E{y i+k .y,} 



(15) 
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Ryx(k)=E{y, + k+ D Xi} 



10 



15 



R yx (k)=E{y i+k .x j } 



M-1 



£ h |- x i+k-l- x i 



M-1 



i=0 



= ^h ( .E{x j+k . l .x j } 



i=0 
M-1 



= ^h,.R xx (k-l) 

t=0 

-Rxx(k)*h k 



(16) 



20 



25 



and similar calculations prove that: 



R yyM=Rxx(k)*(h k *h. k ) 



(17) 



To minimize the mean square error of (13) zero crossings of a partial differential should be found. The target 
impulse response vector B which allows to minimize the mean square error, i.e. the optimal target impulse response 
vector Bopt, in (13) is obtained by solving the equation: 



30 



Substitution of (13) therein results in: 



35 



S(E{e z » 
5(B) 



2R XX .B - 2R Y x-W = 0 



(18) 



(19) 



which implies that the optimal target impulse response vector Bopt can be calculated once the equalizer time domain 
40 vector W is known via the formula: 



Bopt=R^ x . R^.W 



(20) 



45 Herein, R^ 1 x represents the inverse matrix of the above defined autocorrelation matrix Rxx and R"^ x represents the 
transposed matrix of the crosscorrelation matrix Ryx- 

To obtain an expression of the mean square error E which depends only upon the equalizer time domain vector W, 
(20) is substituted in (13). Thus: 



50 
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E{e 2 } = W T .R YY .W + W T R YX .(Rxb T R xx .RiV^x-W (21) 
which is further simplified as follows: 

ITT 1 T T 

(22) 



E{e 2 }= W T (R YY +R YX .(Rxx) T R^ x -2R YX .(Rxx) T R^x) W 



= W T .(R YY -R YX .(Rx 1 x ) T R^ x ).W 
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From Rxx R xx = l ' wherein I represents the unity matrix, it follows that ( R xx) T - R xx = 1 » and that 
( R xx) T = ( R xx) 1 . 

As a result, (22) can be transformed into: 

E{e 2 } = W T (R^ - R Y X'(Rxx)" 1 Rjx)-W (23) 



10 From the definition of the autocorrelation matrix Rxx in (14), it further follows that R^ x = R xx as a resu,t of which: 

E{e 2 } = W T (Ryy - R YX .R x 1 x .Rj x ).W (24 ) 
= W T .O.W 

15 

Herein, O is a channel dependent matrix given by the expression: 

O = R yy " R YX R XX R YX (25) 

20 

When minimizing E{e 2 } by searching the equalizer time domain vector W which minimizes W T . O. W it is obvious that 
the solution: 



25 



W = 



and thus: Bopt = 



(26) 



30 



will be found. However, this solution would imply that no signal could be transmitted over the transmission line. To obtain 
another solution, different from the zero-filter, an energy constraint: 



35 



W.W= 1 



(27) 



has to be fulfilled. 

Thus, the mean square error (24) has to be minimized by taking into account the constraint (27). This is done by solving 
the equations: 



40 



45 



5(E{e 2 }-X.(W T .W-1)) n 

5(E{e 2 }-X.(W T .W.1)) 
6 X 



(28) 



Substituting (24) in (28) and calculating the partial derivative results in: 
50 2.O.W-2.X.W = 0 



(29) 



W'.W-1 =0 

55 The second equation in (29) again represents constraint (27) which has to be fulfilled. The first equation is equivalent to: 

O.W = X.W (30) 
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The equation (30) is an eigenvalue-eigenvector problem. 
Since 

^ min(E{e 2 }) = min(W T .O.W) = min(W T . X.W) 

= min(A,.W T . W) = min(X), 

the optimal equalizer time domain vector Wopt equals the eigenvector of O which corresponds to the minimum eigen- 
value Xmin thereof. 

10 Thus, summarizing the above calculations, the channel dependent matrix O has to be built up, the minimum eigen- 
value Amin and corresponding eigenvector of O have to be calculated. This eigenvector equals the optimized equalizer 
time domain vector Wopt and finally, the optimal target impulse response vector is calculated from Wopt via equation 
(20). Once the matrix O is known, all steps to derive therefrom the equalizer time domain parameters W are thus 
known. 

is Formula (25) shows that the channel dependent matrix O only depends on the transmitted signal X and received 
distorted signal Y. The signal X is transmitted at the initialisation of the system and is well known by the receiver. At 
reception of the distorted signal Y, the receiver is thus enabled to calculate the channel dependent matrix O. The 
dependency between O and the signals X and Y is expressed in (25) via the matrices Rxx, the autocorrelation matrix 
of the transmitted signal X, R Y yi the autocorrelation matrix of the distorted signal Y, and Ryx- the crosscorrelation matrix 

20 of the signals X and Y Notice however that, via Ryx. O also depends on the delay parameter D. This implies that for 
each delay value, a new channel dependent matrix O has to be constructed, each of which giving a new Wopt and Bopt. 

Simplified models and assumptions may allow to simplify the computations necessary to construct the channel 
dependent matrix O. 

In a first simplified implementation of the present invention, the transmitted signal X is supposed to have flat power 
25 spectral density, whilst noise N on the transmission line is supposed to be negligible. Thus, the correlation matrices 
R NN (k) and Rxx( k ) ar © as follows: 

R N N( k )=° 

30 * xx (k)=S x .6(k) 
Therein, 

S x = E(x*). 
5(k)= for k=0, and 

40 6(k)=0 for k * 0. 

The autocorrelation matrix Rxx therefore becomes equal to: 

45 R xx = ^x- 1 ' 

wherein I represents the identity matrix. 

Substituting these expressions for R NN (k) and Rxx(k) in the above expressions (15) for Ryx(k) and RyyM, and substi- 
tuting all the correlation matrices in the channel dependent matrix O (25), it is found that: 



50 



0 = S x H. Txr H TxL .Hj xL (32) 



with: 
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10 



is 



20 



25 



30 



HfxT = 



r h (0) 
r h (T-1) 



- r h (T-1) 

- r h (0) 



(33) 



and: 



T+1 



D+L-T , 



(34) 



In (33): 



M-1 



r h (k)=h k *h. k = ^hi.h^ 

t=0 



(35) 



When measuring in an initial step the line impulse response vector H, the matrix O can be constructed. 

In a second simplified implementation of the present method, the transmission line still is considered to be noise- 
less but the power spectral density of the transmitted signal X is no longer flat. The given power spectral density of the 
signal X has to be taken into account which leads to a new expression for the matrix O. 

In another implementation of the present method, noise N on the transmission line is not negligible but the power 
spectral density of this noise N is known via additional measurements previously executed. In this implementation, a 
new channel dependent matrix O is to be constructed wherein the noise influence is expressed via a term R nn (k) which 
is added to the expression in formula (17). Ryy(k) now is given by: 



RyyC^R^kJ^h^h^+R^k) 



(36) 



35 



40 



45 



50 



Although the above described embodiments exhibit a very low mean square error value and accordingly a negligi- 
ble intersymbol interference outside the guardband, the applicability thereof in ADSL is doubtful. As already indicated 
in the introduction, a concentration of the energy in the unused frequency band has to be avoided. The mean reason 
for his can best be explained by a frequency domain approach. Suppose that H(i) represents the channel transfer func- 
tion in the frequency domain, and E(i) represents an error, equivalent to E, built up in the frequency domain. Thriving 
towards a zero error 



implies that 



E(i)«0 



W(i).Y(i)=B(i).X(i) 



(37) 



(38) 



wherein W(i), Y(i), B(i) and X(i) are frequency domain representations of the equalizer, the distorted signal, the target 
response, and the signal respectively. Equation (38) is equivalent to: 



W(i).^B(i) 



(39) 



55 



and thus to: 



W(i).H(i)=B(i) 



(40) 



since the channel response H(i) is equal to: 
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(41) 



If W(i), H(i) and B(i) are split up into unused band and used band matrices, the above expression (40). because of the 
FDM (Frequency Division Multiplexing) filters used in ADSL is rewritten as: 



w 



r w(i) used 
Lw(i) unu8ed 



'](H(i) unused ] = [ B ^ 



used 



unused 



(42) 



15 



20 



30 



This results in three conditions 

\used 



W(i) used is unconstrained since B(i) used * 0 

W(i) unused . H(i) unused = B(i) unused (43) 
unit energy constraint expressed in the frequency domain 



Although the energy unity norm avoids convergence towards the full zero solution, one notices that the unity norm does 
not guarantee avoidance of a full zero solution in the useful frequency band. Therefore, the equalizer coefficients W 
have to be constrained in order to avoid to concentrate the TEQ energy in the used bandwidth. 
25 In the time domain approach of the above paragraphs, the TEQ energy can be concentrated in the used bandwidth 
by adding an additional degree of freedom. The higher the energy in the unused frequency band, the worse the set of 
equalizer parameters W is chosen. This is expressed by adding the energy in the unused frequency band as an addi- 
tional contribution to the above mean square error E. In this way, one obtains the constrained mean square error CE: 



CE=E+^i.E DC (44) 



Herein, E represents the error contribution due to the distance between the equalized impulse response EIR and the 
target impulse response TIR ( and E DC represents the energy in the unused frequency band or 'don't care' band/The 
35 suppression factor |i is a factor that determines how much the energy in this unused frequency band has to be sup- 
pressed. 

Writing E DC in a matrix notation, one obtains: 
40 E dc=W T . F T DC . F DC .W=W T D. W (45) 

with 



45 



Foc = 



50 



1 e N 



1 e N 



> N 



e 



■j.2.Jr.b4T-1i 
N 



(46) 



Herein, a and b are indexes belonging to the unused or 'don't care* band. 
If (13) and (45) are integrated in (44), following result is found: 



55 



Edc= wT ( r yy + ^- d )+ b T r xx- B " 2 - B T - R YX w 



(47) 



It is noticed that the matrix Ryv in (13) has been replaced by Ryy + fiD. As a consequence, the eigenvector solution 
representing the optimal set of equalizer parameters is found by solving the eigenvalue-eigenvector problem for a new 
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channel dependent matrix O': 

0 , -(R Y y + »i.D)-R YX .R; 1 x .Rjx 

5 

The additional term \x.D herein can be shown to be equivalent to the insertion of a constant virtual noise power spectral 
density component into the unused frequency band. The term 'virtual* is used here to indicate that this is not real meas- 
ured but mathematically injected noise. 

It is remarked here that although the described embodiment in the above paragraphs is a time domain realization, 

10 the present invention also is applicable in frequency domain approaches of the mean square error method. The error 
e k can be rewritten in terms of circular convolutions instead of linear convolutions when using the FFT transform. In this 
way, a frequency domain error criterion is obtained. Such a frequency domain approach will not be described in further 
detail in the present application. However, via the above described time domain example it is clear to a person skilled 
in the art how to adapt e.g. the frequency domain approach of J.S. Chow,J.M. Cioffi and J.A.C. Bingham in their article 

is 'Equalizer Training Algorithms for Multicamer Modulation Systems', Proo. IEEE, p. 761-765, 1993 to be in accordance 
with the present invention. Indeed, such a frequency domain approach can also be adjusted to the FDM requirements 
by modifying the frequency domain mean square error E(i) in a frequency domain constrained mean square error CE(i): 



20 



CE(i)=E(i)+^(i).|W(i) 2 | (49) 



with 



25 _ f // if i belongs to the unused band 

\ 0 if i belongs to the used band 



30 

It is further remarked that, for implementation issues, it is better to avoid injecting noise close to the FDM filter 
boundaries in order to avoid too much attenuation of the first tones of the used band. It is therefore recommended to 
leave a buffer zone without virtually injected noise between the unused and used frequency bands. As an alternative, 
one could use a shaped virtual noise instead of a constant virtual noise. 

35 Still another remark is that intelligent implementations of the present method are taking into account that the calcu- 
lated equalizer coefficients will be rounded. Therefore it is better to reduce the dynamic range and thus to have higher 
gains in the used bandwidth. As a consequence, when different criterions provide equivalent solutions for Wopt and 
Bopt, precision requirements will be taken into account to select one of the solutions. By selecting the solution with least 
dynamic range, the precision performance will be increased. 

40 As already noticed above, the channel dependent matrix O depends on the delay parameter D. Thus, the delay 
value which minimizes the intersymbol interference has to be found to select the best Wopt. The complete theoretical 
search range for D would be from zero to M. However, without significant loss of performance, the search range for D 
can be reduced. To reduce the processing time, one could search for example for D values within a limited interval 
around Dpeak, Dpeak being the index number (between 0 and M) of the peak value of the transmission line impulse 

45 response H. 

Furthermore, it is noted that, to reduce the computation time, the power method can be used to calculate the only 
eigenvector Wopt which has to be found in the present method. This very fast and iterative method is well known but 
extracts only the eigenvector belonging to the largest eigenvalue of a matrix. The above eigenvalue-eigenvector prob- 
lem (30) therefore is rewritten as follows: 

50 

O.W=X.W <=> 0'\W=X" 1 .W (50) 

Thus, an eigenvector of matrix O according to eigenvalue k is also an eigenvector of matrix O" 1 according to the eigen- 
55 value V 1 . The eigenvector Wopt according to the minimum eigenvalue Xmin of O is thus found as the eigenvector 
according to the maximum eigenvalue amin" 1 of the matrix O* 1 . 

While the principles of the invention have been described above in connection with specific apparatus, it is to be 
clearly understood that this description is made only by way of example and not as a limitation on the scope of the 
invention. 
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Claims 



1 . A method for calculating a set of LW equalizer parameters (W), LW being a positive integer number, for an equalizer 
included in a receiver to which a transmitter is coupled via a transmission line and which is able to receive a dis- 
5 torted signal (Y) distorted due to transmission over said transmission line of a signal (X), and which is adapted to 

equalize a transmission line impulse response (LIR) of said transmission line and to thereby generate an equalized 
impulse response (EIR) which approximates a target impulse response (TIR) with a predetermined target impulse 
response length and delayed over a predetermined delay (D), said method comprising the steps of: 

10 a. transmitting said signal (X) and determining said distorted signal (Y) in said receiver; 

b. defining a first integer number of sample periods which represents said target impulse response length; 

c. defining a second integer number of sample periods which represents said predetermined delay (D); 

d. generating a channel dependent matrix (O) from said signal (X), said distorted signal (Y), said target impulse 
response length (L) and said predetermined delay (D), said channel dependent matrix (O) being related to an 

15 error function (E) that has to be minimized; 

e. calculating a minimum eigenvalue (Lmin) of said channel dependent matrix (O); 

f. calculating an eigenvector of said channel dependent matrix (O), said eigenvector corresponding to said min- 
imal eigenvalue (Lmin); 

g. generating said set of LW equalizer parameters (W) from said eigenvector, 

20 

characterized in that said error function (E) comprises a first contribution, representing a difference between said 
equalized impulse response (EIR) and said target impulse response (TIR), and a second contribution representing 
an amount of energy (E DC ) transmitted in unused frequency bands. 

25 2. A method according to claim 1 , characterized in that said channel dependent matrix (O) comprises a first contribu- 
tion, corresponding to said first contribution of said error function (E), said first contribution of said channel depend- 
ent matrix (O) being equal to RYY-RYX.RXX* 1 .RYX T , wherein RYY represents an average autocorrelation matrix of 
a sequence of samples of said distorted signal (Y), said sequence having a length equal to LW, wherein RXX rep- 
resents an average autocorrelation matrix of a sequence of samples of said signal (X), said sequence having a 

30 length equal to said target impulse response length, wherein RYX represents an average correlation matrix of said 
sequence of samples of said signal (X) and said sequence of samples of said distorted signal (Y), and wherein 
RXX" 1 represents an inverse matrix of RXX and RYX T represents a transposed matrix of RYX. 

3. A method according to claim 1 , characterized in that, to generate said channel dependent matrix (O) said transmis- 
35 sion line is supposed to be noiseless. 

4. A method according to claim 1, characterized in that, to generate said channel dependent matrix (O), said signal 
(X) applied to said transmission line by said transmitter is supposed to have a flat power spectral density. 

40 5. A method according to claim 1, characterized in that, to generate said channel dependent matrix (O), said trans- 
mission line is not supposed to be noiseless and that said method further includes a step of measuring noise (N) 
on said transmission line before said channel dependent matrix (O) is generated. 

6. A method according to claim 2, characterized in that RYY(k), which is an element in said autocorrelation matrix 
45 RYY whose row index and column index are different from each other over an integer value k, is equal to 

RXX(k).rh(k) + RNN(k), wherein RXX(k) represents an element in said autocorrelation matrix RXX whose row index 
and column index are different over said integer value k, wherein rh(k) represents a convolution of said transmis- 
sion line impulse response (LIR) and said transmission line impulse response (LIR) delayed over k samples, and 
wherein RNN(k) represents an element in a matrix RNN, whose row index and column index are different over said 
so integer value k, said matrix RNN being an average autocorrelation matrix of noise (N) measured on said transmis- 
sion line. 

7. A method according to claim 1, characterized in that said predetermined delay (D) depends upon amplitude and 
index number of a sample of said transmission line impulse response (LIR) with maximum amplitude. 

55 

8. A method according to claim 1 , characterized in that before executing steps f and g, steps c, d, and e are repeatedly 
executed, said second integer number of sample periods being different in successive executions of step c but lay- 
ing within an interval of acceptable delay values, said channel dependent matrix (O) generated in successive exe- 
cutions of step d as a result thereof being different in such a way that a plurality of minimum eigenvalues similar to 
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said minimum eigenvalue (Lmin) is calculated in successive executions of said step e, and that to perform steps f 
and g, only a said channel dependent matrix with smallest said minimum eigenvalue (Lmin) is selected. 

A device for calculating a set of LW equalizer parameters (W), LW being a positive integer number, for an equalizer 
included in a receiver to which a transmitter is coupled via a transmission line and which is able to receive a dis- 
torted signal (Y) distorted due to transmission over said transmission line of a signal (X), and which is adapted to 
equalize a transmission line impulse response (LIR) of said transmission line and to thereby generate an equalized 
impulse response (EIR) which, when delayed over a predetermined delay (D), approximates a target impulse 
response (TIR) with a predetermined target impulse response length, said device including: 



10 



a. a determining means to determine said distorted signal (Y) received upon transmission of said signal (X); 

b. a first definition means to define a first integer number of sample periods which represents said target 
impulse response length; 

c. a second definition means to define a second integer number of sample periods which represents said pre- 
ys determined delay (D) ; 

d. a first generation means to generate a channel dependent matrix (O) from signal (X), said distorted signal 
00, said target impulse response length and said predetermined delay (D), said channel dependent matrix (O) 
being related to an error function (E) which has to be minimized; 

e. a first calculation means to calculate a minimal eigenvalue (Lmin) of said channel dependent matrix (O); 
20 f. a second calculation means to calculate an eigenvector of said channel dependent matrix (O), said eigen- 
vector corresponding to said minimal eigenvalue (Lmin); 

g. a second generation means to generate said set of LW equalizer parameters (W) from said eigenvector, 

characterized in that said error function (E) comprises a first contribution, representing a difference between said 
25 equalized impulse response (EIR) and said target impulse response (TIR), and a second contribution representing 
an amount of energy (E DC ) transmitted in unused frequency bands. 

10. A demodulator (DEMOD) including between an input (Dl) and an output (DO) thereof, the cascade connection of 
an analog to digital converter (ADC), a time domain equalizer (TEQ), a serial to parallel converter (SFC) a fast fou- 

30 rier transform processing unit (FFT), a frequency domain equalizer (FEQ) and a demapper (DM), said analog to 
digital converter (ADC) being provided to convert a disturbed input signal applied to said demodulator input (Dl) into 
a time domain digital signal, said time domain equalizer (TEQ) being provided to digitally filter said time domain dig- 
ital signal and to thereby generate a time domain equalized signal, said serial to parallel converter (SPC) being pro- 
vided to transform said time domain equalized signal into a time domain parallel sequence, said fast fourier 

35 transform processing unit (FFT) being provided to process said time domain parallel sequence and to thereby gen- 
erate a frequency domain parallel sequence, said frequency domain equalizer being provided to digitally filter said 
frequency domain parallel sequence, and said demapper (DM) being provided to demodulate said frequency 
domain parallel sequence and to generate therefrom an output signal applied to said demodulator output (DO), 
said time domain equalizer (TEQ) including a device for calculating a set of LW equalizer parameters (W), LW being 

40 a positive integer number, for said time domain equalizer (TEQ), said demodulator (DEMOD) being coupled to a 
transmitter via a transmission line and being able to receive a distorted signal (Y) distorted due to transmission over 
said transmission line of a signal (X), and which is adapted to equalize a transmission line impulse response (LIR) 
of said transmission line and to thereby generate an equalized impulse response (EIR) which, when delayed over 
a predetermined delay (D), approximates a target impulse response (TIR) with a predetermined target impulse 

45 response length, said device including: 

a. determining means to determine said distorted signal (Y) received upon transmission of said signal (X); 

b. a first definition means to define a first integer number of sample periods which represents said target 
impulse response length; 

50 c. a second definition means to define a second integer number of sample periods which represents said pre- 

determined delay (D) ; 

d. a first generation means to generate a channel dependent matrix (O) from said signal (X), said distorted sig- 
nal (Y), said target impulse response length and said predetermined delay (D), said channel dependent matrix 
(O) being related to an error function (E) which has to be minimized; 
55 e- a first calculation means to calculate a minimal eigenvalue (Lmin) of said channel dependent matrix (O); 

f. a second calculation means to calculate an eigenvector of said channel dependent matrix (O), said eigen- 
vector corresponding to said minimal eigenvalue (Lmin); 

g. a second generation means to generate said set of LW equalizer parameters (W) from said eigenvector, 
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characterized in that said error function (E) comprises a first contribution, representing a difference between said 
equalized impulse response (EIR) and said target impulse response (TIR), and a second contribution representing 
an amount of energy (E DC ) transmitted in unused frequency bands. 
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